X- 932- 74-133 

PREPRINT 

m 'db. 70 : 

^ ' ' 

GEOS-C ORBIT DETERMINATION 
WITH SATELLITE TO SATELLITE TRACKING 


P. ARGENTIERO 
R. GARZA-ROBLES 


(NASA-T11-X-70659) GEOS-C ORBIT 
OETEBHINATICN WITH SATELLITE TO 
TEACKIHG (NASA) — 21 -»p HC $4.25 

30 


N74-26305 

SATELLITE 
CSCL 22C 

Onclas 

G3/30 40972 


MAY 1974 



Reproduced by 

NATIONAL TECHNICAL 
INFORMATION SERVICE 


US Dopertment of Commerce 
Springfield, VA. 22151 


l 

j 


— GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 



J For information concerning availability 
of this document contact: 

Technical Information Division, Code 250 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 

(Telephone 301-982-4488) 



X-932-74-133 


GEOS-C ORBIT DETERMINATION WITH 
SATELLITE TO SATELLITE TRACKING 


P. Argentiero 
and 

R. Garza-Robles 


May 1974 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 


J 



GEOS-C ORBIT DETERMINATION WITH 
SATELLITE TO SATELLITE TRACKING 


ABSTRACT 

This report studies the feasibility of employing satellite to 
satellite tracking in lieu of ground based tracking to satisfy 
the orbit determination requirements of the GEOS-C mis- 
sion. It is shown that with proper estimation procedures 
it is possible to obtain from S. S. T. data a GEOS-C orbit 
whose altitude error averages about 1 meter. The useful- 
ness of this data type for geopotential recovery is also 
indicated. 
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GEOS-C ORBIT DETERMINATION WITH 
SATELLITE TO SATELLITE TRACKING 


INTRODUCTION 

Satellite to satellite tracking (S. S.T.) by means of a synchronous satellite 
has been recognized as a promising data type for satisfying the stringent 
orbit determination demands of the Earth and Ocean Physics Application 
Program. 2,3, 4 j n report we study the feasibility of employing this data 
type in lieu of ground based tracking to satisfy the orbit determination require- 
ments of the GEOS-C mission. It is shown that with proper estimation pro- 
cedures it is possible to obtain from S. S« T. data a GEOS-C orbit whose altitude 
error averages about 1 meter. This result is significant in itself with respect 
to the GEOS-C mission. It is also significant in that it suggests what the proper 
use of S. S. T. data can accomplish in terms of economical and accurate orbit 
determination of future applications satellites. 


ASSUMPTIONS 

The GEOS-C orbit will be circular with an altitude of approximately 850km and 
an inclination of 65°. The period is 102 minutes. For this study a relay satel- 
lite is assumed at synchronous altitude and at 94° west longitude. The receiving 
station is Rosman. The tracking assumptions and error sources are displayed 
in Table 1. 

With regard to range rate sum data the velocity uncertainty of the relay satellite 
is far more critical than position uncertainty. For this reason a serious mis- 
take is introduced if the relay satellite is treated as an ordinary tracking station 
which happens to be at synchronous altitude. Under such an assumption there is 
no doppler shift due to relative motion between the relay satellite and the ground 
based station and hence the range rate sum data is left unbiased. In actuality the 
relay satellite is in a Keplerian orbit. Consequently an initial position error 
quickly propagates into a velocity error which in turn biases the range rate sum 
data. In this report the velocity errors in the relay satellite are correctly 
propagated into orbital errors of the GEOS-C satellite. 

The results shown in the next section were obtained by means of covariance 
analysis techniques. Hence we have assumed that over the range of expected 
errors, perturbations of orbital estimates are approximately linear functions 
of perturbations of the error sources. With the use of covariance analysis it is 
relatively easy to isolate the effect of any subset of error sources on the orbital 
state estimate at any point along the arc. Such an ability is a valuable asset in 
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Table 1 


Tracking and Error Assumptions for Error Analysis of GEOS-C 
Orbit Determination With S. S. T. Data 


Tracking Assumptions 

Data Type 

Range Rate Sum 

Data Rate 

1/min. 

Arc Length 

5 days 

Data Noise 

0. 1 cm/sec 

Systematic Error Assumptions 

Error Source 

Standard Deviation 

Rosman Survey Error 

10 meters in each component 

Data Bias 

0. 1 cm/sec 

Relay Satellite 

300 meters in each position component 

Epoch State 

3 cm/ sec in each velocity component 

GEOS-C Satellite 

50 meters in each position component 

Epoch State 

5 cm/ sec in each velocity component 

Spherical Harmonic Coefficients 

Full difference between S. A. O. 69 field 

to Degree and Order 8 


and G. E. M. 4 field 


determining which error sources can be ignored and which error sources must 
be simultaneously adjusted along with epoch state. The mathematical develop- 
ment of our covariance analysis techniques is given in the appendix. 


RESULTS 

The simplest estimation of GEOS-C state is obtained when die satellite epoch 
vector is adjusted in a weighted least squares mode and all other parameters 
are fixed at nominal values. The covariance matrix of this estimator was ob- 
tained and propagated along the 5 day arc. The unadjusted error sources and 
associated standard deviations are those shown in Table 1. Figure 1, is a plot 
of the resultant standard deviation of the radial component of the GEOS-C orbit. 
The plot was stable over the entire 5 day arc and hence only the results for the 
first 12 hours are shown. 
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ASSUMPTION 

GEOS-C STATE ESTIMATED. ALL OTHER 
ERROR SOURCES IGNORED. 
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Figure 1. GEOS-C Radial Std. Dev. vs. Time From Epoch 


As shown in the Appendix, covariance analysis techniques permit one to de- 
termine at any point in the arc, the R. S. S. contribution to the standard deviation 
of the estimate of a given state component due to a given unadjusted error source. 
These contributions are called aliasing terms . Table 2 shows the maximum 
aliasing along the arc in the radial, along track, and cross track directions due 
to data bias, Rosman survey error, relay satellite epoch state error, and geo- 
potential. For any given point in the arc and for any given component, the 
aliasing due to geopotential error was obtained by computing the root sum square 
of the corresponding aliasing terms due to errors in individual geopotential 
coefficients up to and including degree and order eight. Similarly the aliasing 
due to relay satellite orbital error is the root sum square of six alias terms 
and the aliasing due to Rosman survey error is the root sum square of three 
alias terms. Table 2 indicates that the dominant error source for GEOS-C 
orbit determination is the uncertainty of the relay satellite state and that data 
bias and Rosman survey error are not significant error sources. 

With regard to GEOS-C orbit determination, attention is focused primarily on 
the altitude estimation requirements imposed by its on-board altimeter. Since 
the instrument is capable of one meter precision it is required to determine the 
GEOS-C altitude with a standard deviation of one meter. From Table 2 it is 
clear that for this to be accomplished the relay satellite state must be simul- 
taneously estimated along with GEOS-C state. Figure 2 shows the resultant 
standard deviation in the estimation of the radial component of GEOS-C position 
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RADIAL STD. DEV. (IN METERS) 


Table 2 


Maximum Aliasing in Meters in Radial, Along Track, and Cross Track 
Components of GEOS-C State Due to Data Bias, Rosman Survey Error, 
Relay Satellite Epoch State Error, and Geopotential Error 



Data Bias 

Rosman 
Survey Error 

Relay Satellite 
Epoch State Error 

Geopotential 

Error 

Radial 

0.025 

0.002 

229 

7.3 

Along Track 

3.300 

0. 008 

1691 

27.3 

Cross Track 

0.210 

0.002 

4379 

17.3 



Figure 2. GEOS-C Radial Std. Dev. vs. Time From Epoch 
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for the first 12 hours after epoch when GEOS-C state is estimated simultaneously 
along with relay satellite state. The average radial uncertainty along the 5 day 
arc was 6.2 meters. 

The object is to satisfy the constraint of one meter radial uncertainty by adding 
the smallest possible number of coefficients to the adjusted parameter set which 
will accomplish the task. Since it is not possible to examine the effects of plac- 
ing each of the possible combinations of coefficients in the adjusted mode, an 
essentially suboptimal procedure for determining an efficient estimation strategy 
must be implemented. Our procedure is recursive. We choose the N + 1st 
coefficient to be adjusted by examining the results of adjusting N coefficients 
along with GEOS-C and relay satellite state and identifying as the next coefficient 
to be adjusted the one which produces the maximum radial aliasing along the 5 
day arc. The recursive process was automated in our covariance analysis 
computer program and was terminated when an average radial error of one meter 
along the 5 day arc was achieved. 

The ordering of coefficients according to their maximum radial aliasing varies 
as terms are added to the adjusted state. This is why the elaborate recursive 
procedure was necessary to identify an efficient estimation strategy. As an ex- 
ample, when geopotential coefficients S(7, 2), S(6,6), C(7, 0), C(7,6), andC(6,5) 
are adjusted along with GEOS-C and relay satellite state the dominant error 
source is the uncertainty of C(7, 2) which provides a maximum radial aliasing of 
2. 9 meters. The uncertainty of geopotential coefficient S(4, 3) is the sixth most 
important error source with a maximum radial aliasing of 1. 9 meters. When 
coefficient C(7, 2) is added to the adjusted state, the uncertainty of coefficient 
S(4,3) becomes the most significant error source with a maximum radial alias- 
ing of 3. 6 meters. Notice also that our recursive procedure was oriented 
toward identifying the most efficient estimation algorithm for minimizing the 
average radial error. If the along track or cross track or the total position 
errors were to be minimized a different set of geopotential coefficients would 
have been identified. 

Table 3 displays the average error in die estimation of radial, along track, and 
cross track components of GEOS-C state as a function of the number of geopo- 
tential coefficients simultaneously adjusted along with GEOS-C and relay satel- 
lite state. Figures 3, 4, and 5 provide the standard deviations of the radial, 
along track, and cross track components of GEOS-C state when the 31 dominant 
geopotential coefficients are adjusted. The 31 coefficients which when adjusted 
yield an average error of one meter in the estimate of the radial component of 
GEOS-C state are given below in the order in which they were identified in the 
recursive procedure: S(7, 2), S(6,6), C(7,0), C(7, 6), C<6,5), C(7, 2), S(4, 3), S{5,5), 
C(4, 4), C(4, 3), S(8,l), C(6, 6), C(8, 7), C(6,4), C(5,5), C(2,2), S<5,4), S(6,4), 
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RADIAL STD. DEV. (IN METERS} 


Table 3 


Average Radial, Along Track, and Cross Track GEOS-C Errors in Meters 
as a Function of Numbers of Estimated Geopotential Coefficients 


Number of Geopotential 
Coefficients Estimated 

Average Radial 
Error (Meters) 

Average Along 
Track Error 
(Meters) 

Average Cross 
Track Error 
(Meters) 

0 

6.2 

43. 7 

63.0 

5 

4.0 

39.6 

60.5 

10 

3.0 

27.2 

42.6 

15 

2.3 

25. 8 

37.3 

20 

1.8 

18.3 

27.0 

25 

1.5 

13.5 

21.3 

31 

1.0 

10.9 

19.5 



Figure 3. Radial Std. Dev. vs. Time From Epoch 
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Figure 4. GEOS-C Along Track Std. Dev. vs. Time From Epoch 


C(7|4), S(6, 5), C(7 , 7), C(5,3), S(4»2), S(7,3), S(7, 7), S(5,3), C(5.1), C(8,4), 
S(6, 2), C(3, 2), S(5, 1). The resultant estimates of these coefficients are rel- 
atively independent and are of good quality. The a priori standard deviations 
of the coefficients were taken to be the difference in representations in the 
S. A. O. 69 field and the G, E. M. 4 field. The resultant standard deviations of 
the estimates of the coefficients averaged about 20% of the a priori standard de- 
viations. This suggests that satellite to satellite tracking is a useful data type 
for geopotential determination. The results also indicate that this estimation 
procedure significantly improves knowledge of relay satellite state. The total 
position error decreased from about 500 meters to about 150 meters. The total 
velocity error decreased from about 5 cm/ sec to about 1 cm/sec. 


CONCLUSIONS 

This paper has shown the feasibility of determining application satellite orbits 
with satellite to satellite tracking. The uncertainty of the relay satellite state 
is the dominant error source. The relay satellite state must be simultaneously 
adjusted along with user satellite state in order to obtain satisfactory results. 
Relay station survey error and data bias are not significant error sources. 
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TIME (IN HOURS) FROM EPOCH 

Figure 5. GEOS-C Cross Track Std. Dev. vs. Time From Epoch 

A covariance analysis was performed to determine the feasibility of determining 
the GEOS-C orbit from 5 days of range rate sum data relayed from a single 
synchronous satellite. It was shown that if relay satellite state and 31 dominant 
geopotential coefficients are simultaneously adjusted along with GEOS-C state 
the radial component of state can be determined with one meter accuracy. The 
estimates of dominant geopotential coefficients as well as relay satellite state 
are also substantially improved. 
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APPENDIX 


COVARIANCE ANALYSIS AS APPLIED TO ORBIT DETERMINATION 


COMPUTING COVARIANCE MATRICES 

Let ^(m) be an m dimensional vector consisting of the differences between the 
correct values of observations of a satellite and nominal values of the observa- 
tions as determined from a nominal orbit. Also let z (n) be an n dimensional 
vector of differences between actual and nominal values of the state of the 
satellite at an epoch and differences between actual and nominal values of 
parameters in the dynamic and measurement models whose associated un- 
certainties may limit our ability to estimate satellite state from the data. The 
sensitivity matrix c (m, n) is defined as that matrix whose element in the ith 
row and the jth column is the partial derivative of y(i) with respect to z(j). A 
first order Taylor series expansion of the functional relationship between y and 
z about the nominal value of z yields 


y - cz (A-l) 


An orbit determination program in processing observations y of y to obtain a 
least square adjustment to z computes a so-called normal matrix defined as 

77 (n, n) * c T wc (A- 2 ) 

where w is a weighting matrix and is usually the inverse of the covariance 
matrix of the observations y of y. Once an orbit determination program com- 
putes and stores the normal matrix, a number of questions can be raised and 
answered at very little cost in terms of computation time. 

The best estimate of the state of the satellite at epoch is obtained by perform- 
ing a least squares adjustment of the state at epoch and all other parameters 
with which are associated significant uncertainties. But frequently this straight- 
forward approach leads to severe core storage requirements. In practice some 
of the parameters in the dynamic and measurement models are estimated along 
with state and others are fixed at their nominal values and left unadjusted in the 
least squares process. In order to determine the consequences of estimating 
some parameters and ignoring others it is useful to compute the covariance 
matrix of such a least squares estimation procedure. 
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Let z be decomposed into two disjoint parameter sets as follows 


Z 


Xi(n,) 


x 2 (n 2 ) 


(A-3) 


where x 1 , is a set of n x , parameters which are to be estimated in a least squares 
process and x 2 is a set of n 2 parameters whose nominal values are left un- 
adjusted by the least squares process but whose uncertainties are to be con- 
sidered in computing the covariance matrix of the resulting estimator. Define 
a matrix A(m, n x ) as a matrix whose element in the ith row and jth column is 
the partial derivative of y(i) with respect to x x (j). Analogously define B(m, n 2 ) 
as the matrix whose element in the ith row and jth column is the partial deriva- 
tive of y<i) with respect to x 2 (j). For future reference notice that the normal 
matrix 7 ? of z as computed and stored by an orbit determination program and 
defined by Equation A - 2 can be written as 


A t wA A t wB 
LB t wA B t wB 


(A-4) 


Assume that there exists a priori estimates of x x and x 2 with properties 


Xj - Xj + aj , E(a.j) “ E(aa T ) - Pj 


x 2 + a 2 , E(a 2 ) _ 0, E(a 2 a 2 ) - P 2 


and assume that the observation vector y or y has properties 


y ~ y + v, E(v) - 0, E(j/v t ) - w 


- 1 


The least squares estimate of x x is obtained as the value of xi which minimizes 
the loss function 


L ( x i) - (y - Ax x - Bx 2 ) T w(y - Axj - Bx 2 ) + (Xj - Xj) 7 Pj 1 ( x i' - Xi) (A-5) 
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The resulting least squares estimator of x t is well known to be 

x t - (A T wA + P- 1 )" 1 [A T w(y - Bx a ') + P^ 1 x/] (A-6) 

Define 


P = [E (x, - x t ) (x t - 


(A-7) 


A series of substitutions reveals that 

Xj - x x = (A t wA + Pj" 1 )' 1 (- A T wBa 2 + A J wu + Pj 1 ^) (A-8) 

Equation 8 yields 

P = (A T wA + Pj -1 ) -1 + (A T wA + Pf 1 ) -1 A T wBP 2 B T wA(A T wA + Pj 1 )” 1 (A-9) 

Notice that the right side of Equation 9 can be computed if one has a priori co- 
variance matrices P 2 and P 2 , and the upper right and upper left portions of the 
normal matrix. To determine the covariance matrix of an estimator which 
estimates some subset of 2 other than , all that is necessary is to permute 
the rows and columns of r/ in the appropriate fashion and proceed as before. 
Thus if one assumes that the normal matrix defined by Equation 2 is precom- 
puted it becomes an easy matter to obtain the resultant covariance matrix when 
any subset of the z parameters are estimated in a least squares sense and the 
rest are ignored. 


THE ALIAS MATRIX 

Assume that all the data has the same variance. Hence 

w = (I^o 2 ) -1 (A- 10) 


where cr 2 is the common variance of each data point. Also assume that the a 
priori estimates of the unadjusted parameters are independent. Under this 
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assumption the covariance matrix P 0 of x 2 ' 


can be written as 


P 


2 



(A- 11) 


where a 2 is the a priori variance of the ith unadjusted parameter. Also define 
a matrix K(nu n 2 ) as 


K - ( A t wA) _ 1 A T wB (A- 12) 

With these assumptions Equation 9 yields the following expression for the ith 
diagonal element of P 


n 2 

P(I. I) = (A- 13) 

j = 0 

where Q is the ith diagonal element of the matrix (A T A)” 1 (this assumes that 
diagonal elements of the matrix Pj 1 are relatively small) and 


fii'j = K(i,j), j>l (A- 14) 

The standard deviation of the ith estimated parameter is given by 


cr. 


1 



(A- 15) 


Define the error sensitivity matrix as 


S = 


= ( ^i.j }• 


= 1 , 2 , 


n, 


j = 0, 1, 


n, (A-16) 
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And finally define the Alias Matrix as 


L = So- 


(A- 17) 


where 


cr 



(A- 18) 


The standard deviation of the ith estimated parameter is seen to be the root sum 
square of the terms in the ith row of the alias matrix. The elements in the first 
column of the alias matrix represent the RSS contribution to the standard devia- 
tion of each estimated parameter due to the data noise. The elements in the jth 
column, j > 2, represent the RSS contribution to the standard deviation of each 
estimated parameter due to the j - 1st unadjusted parameter. 

Possession of the alias matrix reveals much of the probability structure of the 
postulated least squares estimator. With this information one can quickly de- 
termine which error sources are significant with regard to the estimation of a 
given parameter. 

Propagating Covariance Matrices 


Equation 9 provides the covariance matrix of the state x t at some specified 
epoch. In many cases it is important to determine how accurately the state 
can be determined at some time other than epoch. In order to do this cor- 
rectly it is necessary to take into proper account uncertainties in dynamic 
parameters. These parameters may be in an estimated mode or in an un- 
adjusted mode and to incorporate their effect one resorts to state transition 
matrices which presumably have been precomputed by an orbit determination 
program. Let x l (T) be the estimated state at time T. Assume as output from 
an orbit determination program the state transition matrices 


v* (T) 


axj (T) 

dx t 


v 2 (T) 


(T) 

dx 2 


(A- 19) 
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If there are no dynamic parameters in the estimation vector x l5 the matrix 
Vj (T) takes on the particularly simple form, 


VjCT) 


8 0 

0 I 


(A-20) 


where 8 is the six by six matrix defined as the partial derivative matrix of the 
state of the satellite at time T with respect to the state of the satellite at epoch. 
If dynamic parameters are included in the estimated state, the off diagonal 
matrices become non- zero and v x (T) assumes a more complicated form. 

The matrix v 2 (T) is the matrix of partial derivatives of the state x x (T) with 
respect to the unadjusted parameters x 2 . If no dynamic parameters are in 
the unadjusted mode, v 2 (T) is the null matrix. A first order Taylor series 
expansion of the function which describes the time evolution of the state x a (T) 
yields 


*i (T) - v, (T)xj + v 2 (T) x 2 (A-21) 

Substituting x x as obtained from Equation 6 for x x and x 0 ' for xj provides the 
best estimate x x (T), of x x (T) 


*i ( T ) " v i (T) xj + v 2 (T) x 2 (A-22) 

The covariance matrix of Xj (T) is given by 

P (T) = v, (T) P vj (T) + v 2 (T) P 2 v J (T) + Vj (T) E [x x x 2 ] v 2 (T) 


+ v 2 (T) E tx 2 xjl vj (T) (A- 23) 


Equation 23 in conjunction with Equations 6 and 9 yields 


P (T) = vj (T) (A t wA + Pj 1 )" 1 Vj T (T) +[ v x (T) (A T w A + P" 1 )" 1 A T w B 

" v 2 (T)] P^VjCT) (A t w A + P' 1 )' 1 A T wB - v 2 (T)] T (A-24) 
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Finally notice that in much the same fashion that Equation 9 was used to de- 
velop an alias matrix at epoch, Equation 24 can be utilized to develop an alias 
matrix for any time T. 


REMARKS 

If one possesses a functioning orbit determination program it becomes a rela- 
tively easy matter to add covariance analysis capability to the system. A com- 
puter program can be written which assumes as input a normal matrix and state 
transition matrices as generated by the orbit determination program. By 
permuting the rows and columns of the normal matrix and completing the 
matrix operations defined by Equation 9, the covariance matrix of a least 
square process which adjusts any subset of the parameters and ignores the 
rest can be computed. An alias matrix can be obtained and significant error 
sources can be identified. By utilizing the precomputed state transition matri- 
cies, the covariance matrix of the estimate of the state can be propagated from 
epoch to any other time. These operations are very simple and they consume 
little computer time. 

Since the normal matrix and state transition matrices are computed once and 
permanently stored, it is possible to investigate a large number of possible 
estimation strategies. This can be done conveniently and cheaply. For many 
applications such a program is a useful and quickly developed addition to an orbit 
determination system. 
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